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Abstract 

The fractional quantum Hall effect (FQHE) occurs at certain magnetic field 
strengths B*(n) in a two-dimensional electron gas of density n at strong mag- 
netic fields perpendicular to the plane of the electron gas. At these magnetic 
fields strengths, the system is incompressible, i.e., there is a finite cost in 
energy for creating charge density fluctuations in the bulk, while the bound- 
ary of the electron gas has gapless modes of density waves. The bulk energy 
gap arises because of the strong electron-electron interactions. While there 
are very good models for infinite homogeneous systems and for the gapless 
excitations of the boundary of the electron gas, computational methods to 
accurately model finite, inhomogeneous systems with more then about ten 
electrons have not been available until very recently. We will here review an 
ensemble density functional approach to studying the ground state of large 
inhomogeneous spin polarized FQHE systems. 



I. INTRODUCTION 



The fractional quantum Hall effect (FQHE) is manifested in a two-dimensional electron 
gas (2DEG) in a strong magnetic field oriented perpendicular to the plane of the electrons 
[IJ. The effect was discovered as a transport anomaly. This is still the 'Hallmark' of the 
effect, even though there are now a host of other phenomena associated with the effect 
which have been studied experimentally. In a transport measurement it is noted that at 
certain strengths B*(n), which depend on the density n of the 2DEG, current can flow 
without any dissipation. That is, there is no voltage drop along the flow of the current. 
At the same time, the Hall voltage perpendicular to both the direction of the current and 
of the magnetic field is observed to attain a quantized value for a small, but finite, range 
of magnetic field or density, depending on which quantity is varied in the experiment. The 
effect is understood to be the result of an excitation gap in the spectrum of an infinite 2DEG 
at these magnetic fields, so that there is a finite cost in energy to be paid for making density 
fluctuations in the system. This means that the 2DEG is incompressible. In general, the 
magnetic field strengths B*(n) at which the quantum Hall transport anomalies are observed 
are related to the density through the filling factor is = 2n£ 2 B n, with £b — \f^cJ{eB) the 
magnetic length. The quantum Hall effect was first discovered |2| at integer filling factors. 
In this integer quantum Hall effect, the energy gap is nothing but the kinetic energy gap 
hu c = heB/(m*c). Later, the fractional quantum Hall effect was discovered ||] at certain 
rational filling factors of the form v = p/q, with p and q relative primes, and q odd. In the 
FQHE, the excitation gap is a consequence of the strong electron-electron interactions. 

Our understanding of the origin of the FQHE started with Laughlin's seminal paper of 
1983 HJ], which dealt with the simplest fractions v — 1/m, with m an odd integer. At these 
values of is, there are on the average m magnetic flux quanta $o — hc/e per electron. In that 
paper, Laughlin constructed a variational wavefunction for spin-polarized systems in strong 
magnetic fields, strong enough that the splitting hu c between the magnetic subbands, or 
Landau levels, can be taken to be infinite. The wavefunction can then be constructed from 
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single-particle states entirely within the lowest Landau level. Laughlin wrote the variational 
wavef unction as 



where m is an odd integer, and Zj = Xj + iyj is the coordinate of the jth electron in complex 
notation. This wavefunction is an eigenstate of angular momentum. Laughlin went on to 
demonstrate that the system having the wavefunction Eq. ([!]) is an incompressible liquid 
with v = 1/m, m odd, and with an energy gap to excitations, and that the elementary 
excitations are fractionally charged quasi-holes or quasi-particles of charge e* = ±e/m. The 
origin of the energy gap can be understood in the so-called pseudo-potential representation 
of the electron-electron interactions ||. Here, the electron-electron interaction V(ri — r^) 
between electrons i and j is decomposed into strengths in relative angular momentum 
channels £ = 0,1,2,..., of the two electrons. For any realistic interaction V(r< — r^), it turns 
out that Vq > V\ > V% > — Consider the case v — 1/3. In this case, the lowest angular 
momentum pseudo-potential that enters into the m = 3 Laughlin description is V\, the 
interaction energy of two electrons of unit relative angular momentum (in units of H) . (Even 
relative angular momenta are not permissible for spin-polarized electron wavef unctions, since 
they have to be anti-symmetric under interchange of electron coordinates.) The Laughlin 
wavefunction is a very cleverly constructed highly correlated state which completely excludes 
unit relative angular momentum between any two electrons, and is furthermore the only state 
which satisfies this property at v = 1/3. Therefore, any excited state must contain some 
electrons with unit relative angular momentum. The energy gap is due to the cost of this, 
and hence is of order V\. Note that because of the nature of the correlations between the 
electrons, which are contained in the factors [z{ — Zj) m , the Laughlin wavefunction cannot 
be expressed as a single Slater determinant of single-particle states in the lowest Landau 
level. Figure [I] depicts the exchange-correlation energy per particle for infinite, homogeneous 
FQHE systems vs. filling factor. The cusps at filling factors v = 1/3,2/5,3/5, and v = 2/3 
have been included to scale, and these filling factors marked by vertical lines for clarity. 
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Note that these cusps are barely visible on this scale, yet they are responsible for all the 
physics of the FQHE! 

More modern theories of the FQHE are based on the so-called composite fermions. This 
idea was originated by Jain, || who noted that the m = 3 Laughlin wavefunction can be 
written as n ( z % — z jf where = l\(zi — Zj) exp — | \ z k\ 2 is the Slater determinant 
wavefunction of a filled lowest Landau level. Because of constraints on the Hilbert space, 
\&i is an exact eigenstate of an interacting electron system at u = 1. Although need not 
in general be the ground state of the interacting system, it is if the Landau level splitting 
hu) c is much greater than the scale of the Coulomb interaction, e 2 / (6q£b), and if the external 
potential is sufficiently well behaved ||, e.g., it is caused by a uniform positive background 
charge density, and we will assume that these conditions are satisfied. Multiplication by 
the Jastrow factor Yl(zi — Zj) 2 makes the total wavefunction \l/ 3 vanish as the cube of the 
separation of two electrons when they approach each other, rather than linearly. In the 
two-dimensional world of the FQHE, a zero of the wavefunction at a fixed point Zq produces 
a phase factor of 2ir when an electron adiabatically encircles Zq. This is equivalent to adding 
a flux-tube containing a single flux quantum $ & t z . This is purely an Aharonov-Bohm 
phase - no magnetic field is added anywhere except at Zq, where the wavefunction vanishes. 
Therefore, the FQHE at m = 3 can be interpreted as a system of electrons at v — 1, 
but with two flux quanta added to the position of each electron. These composite objects, 
electrons plus an integer number of flux quanta, are called composite fermions. So one can 
say that the fractional quantum Hall effect at v = 1/3 is an integer quantum Hall effect 
[y = 1) of composite fermions. This generation of FQHE states from composite fermions 
was subsequently generalized to all FQHE fractions. One way of studying FQHE systems 
of electrons theoretically is to perform a singular gauge transformation [7], which is a gauge 
transformation in which an odd number of flux quanta is added to the position of each 
electron. The transformed wavefunction is thus a wavefunction for composite fermions. In 
a Lagrangian formulation, a term has to be added to the Lagrangian to ensure that the flux 
tubes indeed are located at the electrons. The resulting term in the Lagrangian is called the 



Chern-Simons term, and is well known from earlier topological field theories. 



II. FINITE SYSTEMS 



We have outlined above how the electron-electron interactions in an infinite, homoge- 
neous system produce the excitation gap. It is important to note that these gaps are only 
for excitations in the bulk of the system. When a system is bounded there must be gapless 
excitations located at the boundaries of the system. The following simple argument, due 
to MacDonald ||, demonstrates this point. Consider a finite system in which the chemical 
potential fi lies in the bulk 'charge gap', i.e., we have to pay the price of the energy gap to 
introduce particles to the bulk of the system. Now imagine that the chemical potential is in- 
creased an infinitesimal amount and consider the resulting change in the current density. 
In the bulk, the current density cannot change since S/j, is infinitesimal and cannot overcome 
the energy gap in the bulk. It follows that if there is a change in the current density as a 
response to 5fi, this must be located at the edges of the system. Current conservation also 
requires that if there is a resulting change in the current along the edge, this change must 
be constant along the edge. We can relate the change in current SI to the change in orbital 
magnetization through 



where A is the total area of the system. This relation is nothing but the equation for the 
magnetic moment of a current loop. But we can write 5M in terms of 5fi using a Maxwell 
relation: 
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where N is the total number of electrons. By combining Eqs. (|2]) and ([J) we arrive at 



Since dn/dB is non-zero, it follows that Sl/Sfj, is non-zero, and we conclude that there 
must be gapless excitations of the system, and that these excitations must localized to the 
edges. Since all experimental systems are finite and inhomogeneous, the low-energy proper- 
ties probed by experiments must be determined by the gapless edge excitations. Advances 
in semiconductor nanofabrication technologies have lead to the possibility of manufacturing 
systems which are extremely inhomogeneous, and in practice dominated by edges. As an 
example, recent experiments have even been performed on tiny quantum dots, with about 
30 electrons on them |]9|,[To[|. In order to accurately understand the experiments and inho- 



mogeneous FQHE systems in general, we must have a way of accurately calculating their 
properties. Certain aspects of inhomogeneous FQHE systems have been studied by different 
techniques. For example, field theories can be constructed to study the low-energy limit of 
the gapless edge excitations fTI||. Composite fermion methods have been used in a Hartree 



approximation to study finite FQHE systems [T^,|T3|. In this approach, the Chern-Simons 
term, arising from the singular gauge transformation, and the electron-electron interaction 
are treated in a self-consistent Hartree approximation. The hope is then that the most im- 
portant aspects of the electron-electron correlations are included in this approximation. Near 
v = 1, at which the Slater determinant is the exact ground state, it makes sense to use 
the Hartree-Fock approximation, and the stability of a quantum dot at v = 1 as a function 
of confining potential has been studied in this approximation [fL4,r5f] . The edge structure has 



also been studied using semiclassical methods [|T(|,|T7]] , in which the electron-electron interac 



tion is included at the Hartree level and it is furthermore assumed that all potentials vary on 



a length scale much larger than lg. Beenakker [16], and Chklovskii, Shklovskii and Glazman 



17fl demonstrated that the edge of an integer quantum Hall system, in which the correlation 
energies between the electrons can be ignored, consists of a sequence of compressible and 
incompressible strips. Imagine going from the bulk of an integer quantum Hall effect, with 
v filled Landau levels in the bulk (so that the chemical potential lies above the energy of 
these v Landau levels), toward an edge. There is an external potential confining the system, 
and this potential rises toward the edge, causing the Landau levels to bend upward near the 



edge. The compressible strips occur where the chemical potential crosses a Landau level. 
There are then both empty and occupied single-particle states available, and the electron 
gas can screen the external potential perfectly. Eventually this Landau level bends upward, 
rising above the chemical potential. There are then no more empty single particle states, and 
changing the electron density would involve a cost of energy of the order of hu c , the spacing 
between the Landau levels. Thus, an incompressible strip forms. Here, the electrons cannot 
screen the external potential. Further out on the edge, the chemical potential crosses the 
next lower Landau level, and a new compressible strip forms. Thus, whenever the chemical 
potential lies between two Landau levels, the electron gas is incompressible, and the elec- 
tron density is constant, while the total potential varies. On the other hand, whenever the 
chemical potential crosses a Landau level, the electron gas can screen perfectly the external 
potential, and the electron density varies while the total potential remains constant. The 
width of the incompressible strips is then determined by the length over which the confining 
potential varies an amount equal to the energy gap hu c . The origin of the compressible and 
incompressible strips are the energy gaps, which are the kinetic energy gaps hu c in the case 
of the integer quantum Hall effect. But it is easy to heuristically generalize the argument to 
include the energy gaps causing the FQHE |I7| . The conclusion is then that there should be 



compressible and incompressible strips, with the density of the incompressible strips fixed 
at the density of an FQHE fraction. The width of each incompressible strip is then fixed by 
the length over which the confining potential varies an amount equal to the energy gap of 
the FQHE fraction corresponding to the density of that strip. 

Finite, inhomogeneous systems have also been studied by direct numerical diagonaliza- 



tions |18[ . At the present, numerical diagonalizations are limited to systems with of the 
order of 10 electrons. 

It is highly desirable to have a computational approach which accurately includes 
electron-electron correlations and can handle inhomogeneous systems with on the order 
of 10 2 -10 3 electrons. One such approach which is in principle valid for any interacting elec- 
tron system is the density functional theory (DFT) [|l9H 2l| . There have been some attempts 



to apply density functional theory to the FQHE. Ferconi and Vignale [^] applied current 
density functional theory |23| to small, parabolically confined quantum Hall systems and 
showed that the current density functional theory gave good results for the ground state 
energy and spin polarization near v — 1. However, the energy gaps due to correlation effects 
were not included in that calculation. Ferconi, Geller and Vignale EM also recently studied 
FQHE systems within the spirit of the DFT using an extended Thomas-Fermi approxima- 
tion at low, but non-zero, temperatures. In this, the kinetic energy was treated as a local 
functional, as in the standard Thomas-Fermi approximation, while the exchange-correlation 
energy was included in a local density approximation (LDA). This extended Thomas- Fermi 
approximation is valid in the limit of very slowly varying confining potential. Ferconi, Geller 
and Vignale focused on the incompressible and compressible strips at an edge of an FQHE 
system, and obtained results in agreement with the predictions by Chklovskii, Shklovskii, 
and Glazman [17] . 



We have developed for the fractional quantum Hall effect an ensemble DFT scheme 
within the local density approximation, and have applied it to spin-polarized circularly 
symmetric quantum dots . In our approach, the kinetic energy is treated exactly, and the 
density represented by Kohn-Sham orbitals. The results are in good agreement with results 
obtained by semiclassical |16|JT7| , |2^| , Hartree-Fock |l^JT5|] (for cases where the correlations 
do not play a major role), and exact diagonalization methods ||26|| . Our calculations show 
that the exchange and correlation effects of the FQHE are very well represented by the LDA 
and that our approach provides a computational scheme to model large inhomogeneous 
FQHE systems. We note that there exist previous formal DFTs for strongly correlated 
systems, in particular for high-temperature superconductors |p7H , and DFT calculations of 
high-T c materials |28| and transition-metal oxides |2£|. However, ours are, to the best of 
our knowledge, the first practical LDA-DFT calculations of a strongly correlated system in 
strong magnetic fields, and demonstrate the usefulness of the LDA-DFT in studying large 
inhomogeneous FQHE systems. 



S 



III. ENSEMBLE DENSITY FUNCTIONAL THEORY APPROACH 



In typical DFT calculations of systems of N e \ electrons, the standard Kohn-Sham (KS) 



scheme [30] is implemented, in which the particle density n(v) is expressed in terms of a 
Slater determinant of iV > N e \ KS orbitals, ip a {r). These obey an effective single-particle 
Schrodinger equation H eS ip a = e a ip a , which is solved self-consistently by occupying the 
iV e i KS orbitals with the lowest eigenvalues e a , and iterating. This scheme works well in 
practice for systems for which the true electron density can be represented by a single Slater 
determinant of single-particle wavef unctions. However, when the KS orbitals are degenerate 
at the Fermi energy (which we identify with the largest e a of the occupied orbitals) there 
is an ambiguity in how to occupy these degenerate orbitals. There exists an extension of 
DFT which is formally able to deal with this situation. This extension is called ensemble 



DFT poy2l| , and in it, the density of the system is represented by an ensemble of Slater 
determinants of KS orbitals. However, while it can be shown using ensemble DFT that 
such a representation of the density is rigorous, it cannot be shown how the degenerate KS 
orbitals at the Fermi energy should be occupied, i.e., there has not been available a practical 
computational scheme for ensemble density functional theory. 

We argued above that the Laughlin wavefunction cannot be represented as a Slater 
determinant of single-particle wavefunctions. Therefore, one may suspect that the density 
of a general FQHE system cannot be represented by a single Slater determinant. We will 
now argue that this is indeed the case. Consider a FQHE system in the xy-plane with 
the magnetic field along the z-axis. A circularly symmetric external potential V^ x t(r) = 
K xt (r) (due, e.g., to a uniform positive background charge density) confines the systems 
such that the density is fixed with a local filling of v = 1/3 up to an edge at tq (tq ^> Ib) 
where the density falls to zero within a distance of order £b- That such systems exist 
is well demonstrated by the excellent agreement between the Laughlin wavefunction and 



experiments, and by many numerical calculations [^T|,^]. Due to the circular symmetry 
we can label single-particle orbitals by angular momentum m, and by a Landau level index 
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n > 0. The orbitals ift m ,n{ T ) are centered on circles of radii r m w \/2m£B with Gaussian 
fall-off s for r <C r TO and r ^> r m . The single-particle orbitals with n = are then in the 
bulk all degenerate, and the degeneracy is not lifted by electron-electron interactions since 
the system is homogeneous in the bulk. Only the orbital ?/>o,o is non-zero at the origin - all 
others vanish at r = 0. In order to obtain a constant density at v — 1/3 even at the center 
of the system, all single-particle orbitals in the bulk with n = must have occupancies 
1/3. If the Fermi energy lay above the energies of the bulk orbitals, they would all be 
filled and one would have v — 1. Therefore, to get occupancies 1/3 the Fermi energy 
must lie at the degenerate energy e m o of these orbitals. Thus, in applying DFT to the 
FQHE we can expect a huge degeneracy of KS orbitals at the Fermi energy, and they must 
all have fractional occupancies. Consequently, the particle density cannot be expressed in 
terms of a single Slater determinant. Instead, the density has to be constructed from an 
ensemble of Slater determinants, i.e., the orbitals at the Fermi energy are assigned fractional 
occupation numbers, just as the Laughlin wavefunction for v = 1/3 is not a single Slater 
determinant, but a highly correlated state with average occupancies of 1/3 of single-particle 
states. Therefore, one might expect from the outset that one has to use ensemble density 
functional theory - the standard Kohn-Sham scheme may not converge. 

Although ensemble DFT has been developed formally, there are in practice few examples 
of applications and calculations using ensemble DFT for ground state calculations. A sig- 
nificant aspect of our work is that we have developed an ensemble scheme which is practical 
and useful for the study of the FQHE. In ensemble DFT, any physical density n(r) can 
be represented by n(r) = Ximn fmn\ipmn(j)\ 2 , where f mn are occupation numbers satisfying 
< fmn < 1, and the orbitals ip mn satisfy the equation 

where V x A(r) = B(r). In equation ([5]), Vk(r) is the Hartree interaction of the 2D elec- 
trons, and, as usual, V xc (r, B) is the exchange-correlation potential, defined as a functional 
derivative of the exchange-correlation energy E xc [n(r), B] of the system with respect to den- 
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sity: V xc (r, B) = SE *f™ff' B ^ ■ (We will hereafter not explicitly indicate the parametric 
dependence of V xc and E xc on B.) For the case of the FQHE, we know that the exchange- 
correlation potential will be crucial, as it contains all the effects of the electron correlations 
which cause the FQHE in the first place, and a major part of the DFT application is to come 
up with an accurate model of E xc and so of V xc . Leaving this question aside for a moment, 
and assuming that we have succeeded in doing so, the practical question is then how to de- 
termine the KS orbitals and their occupancies in the presence of degeneracies. We devised 
an empirical scheme, which means that after a lot of trial and error we made some educated 
guesses that work. Our scheme produces a set of occupancies for the KS orbitals which 
satisfy some minimum requirements, namely (a) the scheme converges to physical densities 
(to the best of our knowledge) for FQHE systems, (b) it reproduces finite temperature DFT 
distributions at finite temperatures, and (c) it reproduces the standard Kohn-Sham scheme 
for systems whose densities can be represented by a single Slater determinant. 

In our scheme, we start with input occupancies and single-particle orbitals and iterate the 
system N eq times using the KS scheme. The number N cq is chosen large enough (about 20- 
30 in practical calculations) that the density is close to the final density after N eq iterations. 
If the density of the system could be represented by a single Slater determinant of the KS 
orbitals, we would now essentially be done. However, in this system there are now in general 
many degenerate or near-degenerate orbitals at the Fermi energy. After each iteration, 
the Kohn-Sham scheme chooses to occupy the iV el orbitals with the lowest eigenvalues, 
corresponding to making a distinct Slater determinant of these orbitals. But there will be 
small fluctuations in the density between each iteration, which cause a different subset of 
these (near) degenerate orbitals to be occupied after each iteration. This corresponds to 
constructing different Slater determinants after each iteration, and the occupation numbers 
f mn of these orbital are zero or unity more or less at random after each iteration. This 
means that the computations will never converge. However, the average occupancies, i.e., 
the occupancies averaged over many iterations, become well defined and approach a definite 
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value, e.g., 1/3 for orbitals localized in a region where the local filling factor is close to v — 
1/3. Therefore, we use these average occupancies to construct an ensemble by accumulating 
running average occupancies (f mn ) after the initial N eq iterations 

I N it 

(fmn) = TTT T ^ fmn,i, (6) 

^V it - iV cq j i=Ncq + 1 

where f m n,i is the occupation number (0 or 1) of orbital ip mn after the ith iteration, and 
use these to calculate densities. Thus, our algorithm essentially picks a different (near) 
degenerate Slater determinant after each iteration, and these determinants are all weighted 
equally in the ensemble. It is clear that this scheme reduces to the KS scheme for which the 
density can be represented by a single Slater determinant of KS orbitals (for which the KS 
scheme picks only the one Slater determinant which gives the ground state density) for N eq 
large enough. We have numerically verified that a finite-temperature version of our scheme 
converges to a thermal ensemble at finite temperatures down to temperatures of the order 
of lO~ 3 huj c /kB. We have also performed some Monte Carlo simulations about the ensemble 
obtained by our scheme. In these simulations, we used a Metropolis algorithm to randomly 
change the occupation numbers about our converged solution, keeping the chemical potential 
fixed. The free energy of the new set of occupation numbers was calculated self-consistently. 
If the free energy decreased, this set was kept, and if the free energy increased, the set was 
kept if a random number was smaller than exp [—AF/k B T*j, where AF is the change in free 
energy, and T* a fictitious temperature. The results were that to within numerical accuracy 
our ensembl DFT scheme gives the lowest free energy. As a condition for convergence, we 
typically demanded that the difference between the input and output ensemble densities, 
ni n (r) and n out (r), of one iteration should satisfy 

— \n in (r)-n out (r)\rdr<10- 3 . (7) 

N e \ JO 

Practical density functional theory calculations hinge on the availability of good approx- 
imations for the exchange-correlation potential V xc , which enters in the effective Schrodinger 
equation for the KS orbitals. The simplest, and probably the most commonly used, approx- 
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imation is the local density approximation (LDA). In this approximation, the exchange- 
correlation energy is assumed to be a local function of density, so that the total exchange- 
correlation energy consists of contributions from the local density of the system. Thus, in this 
approximation one writes E xc /N = J dre xc (z/)n(r), where e xc (z/) is the exchange-correlation 
energy per particle in a homogeneous system of constant density n = v/{2n^ B ) and filling 
factor v. In other words, in the LDA one assumes that the system is locally homogeneous, 
i.e., the system can locally be approximated to have the energy per particle of an infinite, 
homogeneous system of the local density. This approximation obviously makes sense if the 
density of the system varies on a very long length scale, while it could be questionable for 
systems in which the density varies on some microscopic length scale. However, experience 
has shown that the LDA often works surprisingly well, even for systems in which the electron 
density is strongly inhomogeneous fT9|. In fact, the first application of LDA-DFT was to 



calculate the work function of simple metals [32], so these were systems which were termi- 



nated with densities varying on the scale of a Bohr radius! Nevertheless, the LDA-DFT gave 
quite good results, vastly superior to those of the Hartree- or Hartree-Fock approximation. 
In the FQHE, the length scale of exchange-correlation interactions and density fluctuations 
is given by the magnetic length is due to the Gaussian fall-off of any single-particle basis in 
which the interacting Hamiltonian is expanded. The densities are relatively smooth on this 
length scale, which gives us additional hope that the LDA will work well for the FQHE, too. 
In addition, the cusps in the exchange-correlation energy will suppress density fluctuations, 
so in this sense one can actually expect the basic physics of the FQHE to make the LDA a 
good approximation. 

In conventional LDA-DFT calculation, the exchange-correlation energy e xc is obtained 
by interpolating between the exchange-correlation energies per particle of systems with van- 
ishing and infinite densities, respectively, for which exact results are known. Analogously, 
following Rasolt and Perrot ||33|| , we obtain our exchange-correlation energy by interpolating 
between two limits for which the result is known very accurately. In our case, the two limits 
are B — > oo, and B — > 0, respectively, and we stitch them together using a Pade approx- 
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imant [ 33 1 . Thus, we write for the exchange-correlation energy per particle of a uniform 



electron gas in a constant magnetic field 

M^J ' ^ } 

Here, ej c c is the zero-magnetic field result for a 2DEG obtained by Tanatar and Ceperley 
3~4fl . The term e^ HE (z/) is the B — > oo limit, which is the exchange-correlation energy of 



the FQHE in a system for which only single-particle states in the lowest Landau level are 
occupied. This contribution consists of two terms. The first one is a smooth interpolation 
formula e^ M (z/) due to Levesque, Weiss and MacDonald |35|] between ground state energies 



at some rational fillings. The second one, e^ c (V), is all- important for the study of the 
FQHE. This term contains the cusps in the ground state energy which cause the FQHE. 
Here we have used a simple model which captures the essential physics. We model e^ c (u) 
by constructing it to be zero at values of v = p/q which display the FQHE. Near v = p/q, 
^ c (v) is linear and has at v = p/q a discontinuity in the slope related to the chemical 
potential gap A/i = g(|A p | + \A h \). Here A Pth are the quasiparticle (hole) creation energies 
which can be obtained from the literature at fractions v = p/q. Farther away from 

v = p/q, e^ c {y) decays to zero. Finally, in the LDA V xc (r) is obtained from e xc (z/) as 
V xc (r) = 9[t/e , xcM1 at constant B. In our calculations, we restrict ourselves to include 

only the cusps at v — 1/3,2/5,3/5 and v = 2/3, which are the strongest fractions. These 
are some the fractions of the form v = ^ 2 p±i) g enera t e d by the so-called Vj-model, in which 
only the pseudo-potential V\ is included. 

A technical difficulty arises in the LDA: the discontinuities in V xc (r) in the LDA give rise 
to a numerical instability. The reason is that an arbitrarily small fluctuation in charge density 
close to an FQHE fraction gives rise to a finite change in energy. Imagine that the local 
filling factor z/(r) in some neighborhood of a point r is very close to, but less than, say, 1/3 
after one iteration. In this neighborhood, the local exchange-correlation potential will then 
form a potential well with sharp barriers at the points around r where z/(r) = 1/3. During 
the next iteration, charge will then be poured into this well. As a result, the local filling 
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factor will after this iteration exceed 1/3, and in this neighborhood V xc now forms a potential 
barrier of finite height. So in the next iteration, charge is removed from this neighborhood, 
and so on. We can see that this leads to serious convergence problems. To overcome this, 
we made the compressibility of the system finite, but very small, corresponding to a finite, 
but very large, curvature instead of a point-like cusp in e xc at the FQHE fractions. In 
other words, instead of having a step-like discontinuity A/x in the chemical potential, it rises 
smoothly an amount A/i over an interval 7 in the filling factor. What we found worked very 
well in practice was to have the discontinuity in chemical potential occur over an interval of 
filling factor 7 of magnitude 10~ 3 . This corresponds to a sound velocity of about 10 6 m/s 
in the electron gas, which is three orders of magnitude larger than the Fermi velocity of a 
2D electron gas at densities typical for the FQHE. In general, the finite compressibility does 
not lead to any spurious physical effects so long as the energy of density fluctuations on a 
size of the order of the systems size is larger than any other relevant energy in the problem. 
The only noticeable effect is that incompressible plateaus, at which the density would be 
perfectly constant were the compressibility zero, will have density fluctuations on a scale of 
7. Figure |2] depicts V xc used in our calculations as a function of filling factor. 



IV. APPLICATIONS TO QUANTUM DOTS 

We have self-consistently solved the KS equations Eqs. (||) for a spin-polarized quantum 
dot in a parabolic external potential, V ext (r) = ^m*fl 2 r 2 , by expanding the KS orbitals 
ipmn( r ) = e irn( ^(p mn (r) in the eigenstates of H = (p + ~A(r)) . We use the cylindrical 
gauge, A(r) = and include the four lowest Landau levels (n = 0, ... ,3). We chose 

the static dielectric constant eo = 13.6, appropriate for GaAs, and a confining potential of 
strength |J HQ = 1.6 meV. 

The use of our LDA-DFT scheme is illustrated by a study of the edge reconstruction of 
the quantum dot as a function of magnetic field strength. As is known from Hartree-Fock 
and exact diagonalizations P||I1^P^,|T2"|.P"3] , for strong confinement the quantum dot forms 
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a maximum density droplet in which the density is uniform at v — 1 in the interior, and 
falls off rapidly to zero at r ~ V2N£b = r . As the magnetic field strength increases, 
a "lump" of density breaks off, leaving a "hole" or deficit at about r = r$. This effect 
is due to the short-ranged attractive exchange interaction: it is energetically favorable to 
have a lump of density break off so that the system can take advantage of the exchange 
energy in the lump. As B is further increased, the correlations will cause incompressible 
strips with densities v = p/q to appear |T6| , p!7|j38| , |2^| on the edges, and incompressible 



droplets to form in the bulk at densities u = p/q. Figure || depicts various stages of edge 
reconstruction obtained by us as the magnetic field strength is increased. The value of 
B for which the exchange lump appears compares very well with the value found by De 
Chamon and Wen |15|] in Hartree-Fock and numerical diagonalizations. At higher fields 



still, incompressible strips appear at the edges, and incompressible droplets are formed 
in the bulk. In figure (| we show occupancies for the KS orbitals for a finite-temperature 
calculation with N e \ = 40, B = 4.45 T, and T = 0.003 e 2 / (e ^B^s)- The diamonds depict the 
converged ensemble occupancies for the KS single-particle states using our ensemble scheme. 
At this finite temperature, we calculated the thermal occupancies of the KS orbitals after 
each iteration, and these 'instantaneous' thermal occupancies were then averaged using our 
ensemble scheme. The temperature was sufficiently low that the 'instantaneous' thermal 
occupancies were essentially or 1 before convergence. The continuous curve shows the 
'instantaneous' thermal occupancies obtained after a particular iteration after convergence 
has been achieved. This figures then clearly shows that our ensemble occupancies (in this 
case at a low, but nonzero, temperature) converged to the thermal occupancies. Note that 
this particular temperature is so low that no standard finite-temperature scheme could be 
used to achieve convergence. 

Figure |5] depicts the eigenvalues of the KS orbitals for iV el = 40, and B = 4.45 T. The 
dashed line indicates the chemical potential of the system. This figure then shows that 
all KS orbitals in the bulk are in fact degenerate. It may at first seem paradoxical that 
the eigenvalues are degenerate on an incompressible strip, since, according to the picture 
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by Chklovskii, Shklovskii, and Glazman ||17|| , on such a strip the density is constant, while 
the total potential varies (since the electrons cannot screen the external potential). If the 
total potential varies, then ought not the the eigenvalues of the KS orbitals localized on 
that strip vary, too, since these then in general are subjected to different potential energies? 
The problem with this argument as applied to DFT is that it ignores the effect of the 
exchange-correlation potential. As the external and Hartree potentials vary across the strip, 
the exchange-correlation potential varies across its discontinuity so as to completely screen 
out the external and Hartree potentials. The discontinuity in V xc does not mean that this 
potential is fixed at the lower limit of its discontinuity while the density is fixed at an 
incompressible strip. What it does mean, is that V xc is free to achieve any value across its 
discontinuity so as to completely screen out the external and Hartree potentials. In this 
way, it is perhaps better to think of incompressibility as the limit of a finite compressibility 
approaching zero. A strip can then remain incompressible with constant density so long as 
V xc can screen the external and Hartree potentials, so the width of the incompressible strip 
is given by the distance over which the external plus Hartree potentials varies an amount 
given by the energy gap associated with the density at that strip. Also, all bulk KS states 
are degenerate at the chemical potential. When a single particle is added, the chemical 
potential simply increases a small amount, and all KS orbitals are again degenerate at the 
chemical potential. We also would like to emphasize that incompressible regions that appear 
in these calculations are not due to the presence of a uniform positive background density 
which tends to fix the bulk density at the value of the background density. 

There is also another edge effect caused by correlations. For particular, stiff confining 
potentials, so-called composite edges 0,^] can appear. These can be thought of as particle- 
hole conjugates of uniform incompressible droplets. Consider a droplet with a bulk density 
corresponding to v = 1/3, falling off to zero at the edge. A incompressible droplet with a 
bulk density of v = 2/3 is obtained by particle- hole conjugation. However, at the edge, the 
density will first rise to v = 1 (since the density of the v = 1/3 droplet drops to zero), and 
then eventually drop to zero. Note that this argument is based on particle-hole conjugation, 
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which is an exact symmetry of the lowest Landau level |40| , and it is unclear if composite 



edges exist in real systems, which do not strictly obey particle-hole symmetry. 

Figure |6] depicts the particle density (inset) for a system where the confining potential 
is supplied by a uniform positive background charge density n + = 2/(6ii£ 2 B ) (so that the 
corresponding filling factor is u + = 2/3) for r < r , and falling linearly to zero over a distance 
a for a > tq, where ro is fixed by charge neutrality. Thus, the parameter a is a convenient 
parameter with which one can control the stiffness of the confining potential . From this 
figure, we see that for a = 0, the system forms a composite edge, even though our system 
does not obey particle-hole symmetry. We therefore conclude that such edges can exist in 
real systems. As a is increased, the edge undergoes an instability and reconstruction, and 
eventually forms incompressible strips. 

V. CONCLUSION AND SUMMARY 

In conclusion, we have showed that ensemble density functional theory can be applied 
to the FQHE. This opens the door to doing realistic calculations for large systems. We 
believe that our results are also significant in that they are the first LDA-DFT calculations 
of a strongly correlated system in a strong magnetic field, and they are (to the best of our 
knowledge) the first practical ensemble DFT calculations. There are, however, still many 
issues that need to be resolved, and new directions to go. For example, our calculations 
were of a spin-polarized system. As is well known ||, the spin degree of freedom is very 
important, even for magnetic fields of the order of 10 T. The reason is that the effective 
g-f actor in GaAs is very small, so small that the ratio between the Zeeman splitting and the 
cyclotron energy is about 2%. This leads to the possibility of FQHE ground states which 
are not spin-polarized. It also gives rise to so-called charge-spin-texture excitations near 
v = l/m, m — 1,3, .... These excitations are lower in energy than a simple singlet particle- 
hole pair, and believed to be responsible for the observed |^TJ rapid destruction of the spin 



polarization near v = 1. Correctly including the spin degree of freedom to account for this 
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involves a DFT for Heisenberg spins. We are presently, together with J. Kinaret (Chalmers 
University of Technology) developing such a theory. 
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FIGURES 

FIG. 1. The ground state energy e X c per particle of an infinite, homogeneous, spin-polarized 
FQHE system is depicted as a function of filling factor. The cusps at v = 1/3, 2/5, 3/5, and v = 2/3 
are included (these filling factors are indicated by vertical lines for clarity). 

FIG. 2. Exchange-correlation potential V xc as function of filling factor in units of e 2 /(eo^s) for 
< v < 1. The increase in V xc at an FQHE filling occurs over a range of filling factor of 0.004. 

FIG. 3. Edge reconstruction of a quantum dot as the magnetic field strength is increased. 
Plotted here is the local filling factor v(r) for a parabolic quantum dot with HQ = 1.6 meV and 40 
electrons. For magnetic field strengths B<2.5T the dot forms a maximum density droplet, and 
for B m 3.0 T, an exchange hole is formed. For stronger magnetic fields, incompressible regions 
form, separated by compressible strips. 

FIG. 4. Ensemble (diamonds) and 'instantaneous' (solid line) thermal occupancies for 7V e i = 40, 
B = 4.45T, and T = 0.003 e 2 /(e^ B k B ) after convergence. 

FIG. 5. Eigenvalues of the lowest-Landau level Kohn-Sham orbitals for N c \ = 40 and B = 4.45 
T as a function of angular momentum quantum number. The chemical potential is indicated by 
the dashed line. 

FIG. 6. Local filling factor v(r) as a function of r (in units of Ib) for a system of 45 electrons 
in a magnetic field of B = 5.0 T. The confining potential is due to a positive background charge 
density at v + = 2/3 in the bulk, and falling linearly to zero within a distance a near the edge. 
For a stiff edge (a = 0), the system forms a composite edge with v(r) rising towards unity near 
the edge. As a increases, the edge becomes softer and undergoes a reconstruction to a sequence of 
incompressible strips. 
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